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The spin structure of the pion 
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We present the first calculation of the transverse spin structure of the pion in lattice QCD. 
Our simulations are based on two flavors of non-perturbatively improved Wilson fermions, with 
pion masses as low as 400 MeV in volumes up to (2.1 fm)^ and lattice spacings below 0.1 fm. We 
find a characteristic asymmetry in the spatial distribution of transversely polarized quarks. This 
asymmetry is very similar in magnitude to the analogous asymmetry we previously obtained for 
quarks in the nucleon. Our results support the hypothesis that all Boer- Mulders functions are alike. 
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Introduction. — Since their discovery in the late 1940s, 
pions have played a central role in nuclear and particle 
physics. As pseudo-Goldstone bosons of spontaneously 
broken chiral symmetry they are at the core of the low- 
energy sector of quantum chromodynamics (QCD). Since 
the pion has spin zero, its longitudinal spin structure in 
terms of quark and gluon degrees of freedom is trivial. 
Pion matrix elements of quark and gluon helicity opera- 
tors vanish due to parity invariance, (7r(P')|I]'^|7r(P)) =0, 
where, e.g., for quarks T,'^ = qj^^5q. An instructive quan- 
tity describing the spin structure of hadrons is the prob- 
ability density p{x,b±) of quarks in impact parameter 
space [ij , illustrated in Fig. |T] 
Here x is the longitudinal mo- 
mentum fraction carried by 
the quark, and the impact pa- 
rameter b± gives the distance 
between the quark and the 
center of momentum of the 
hadron in the plane transverse 
to its motion. Because of 
parity invariance, the density 
p{x,bj^,X) of quarks with he- 
licity A in a pion is determined 
by the unpolarized density, 
2p{x,b±,X) = p{x,b±). The 
latter is given by p{x,b±) — 
H"^ {x, £^=0,b^) in terms of a 

b± dependent generalized parton distribution (GPD) at 
zero skewness The lattice QCD calculations discussed 
below give access to x-moments of quark spin densi- 
ties, which we have investigated in 3] for quarks with 
transverse spin s± in a nucleon with transverse spin Sj_ . 
The corresponding expression p(x,b±, s±) for polarized 
quarks in the pion is obtained by setting S± = in the 
nucleon densities of [1, The result is much simpler 



FIG. 1: Illustration of 
the quark distribution in 
a TT^ in impact parameter 
space. 



but still contains a dipole term cx s'j^e'-' b^, which leads 
to a dependence on the direction of b± for fixed s±, 



P"(6i,s±) 



dxx"" ^p{x,bj_,s±^ 



,(1) 



where B 



TnO 



9f,2^ _Bj^,jQ . The b± dependent vector and 
tensor generalized form factors (GFFs) of the pion, A^q 
and i?J„o' ^''^ moments of the GPDs: 



/. 



dxx''-^H''{x,^^0,bl} 



dx ; 



'E^ix,^=0,bl) - B^^bl). 



(2) 



To this day, next to nothing is known about the signs 
and sizes of the B^^q. Since these GFFs determine the 
dipole-like distortion of the quark density in the trans- 
verse plane, non- vanishing -B^'no would imply a surprising 
non-trivial transverse spin structure of the pion. A com- 
putation of the B^^Q from first principles in lattice QCD 
therefore provides crucial insight into the pion structure. 

Lattice QCD calculations give access to GFFs F{t) — 
A^o{t) , B^no(^) momentum space, which are related 
to the impact parameter dependent GFFs P(fo5_) = 
^no(^i)' ^Tno(^l) ^7 a Fouricr transformation 



d^A±e 



-ib I -Aj 



F(< = -Ai), (3) 



where is the transverse momentum transfer. The 
momentum-space GFFs B^^Q(t) parameterize pion ma- 
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FIG. 2: Lattice results at /3 = 5.29 and m-n ~ 600 MeV for 
the first two generalized form factors B^'J^Q{t) for up-quarks 
in the n'^ . The shaded bands show p-pole parameterizations. 

trix elements of local tensor quark operators, 



ji-i 



^ A^i • • • A^-P^-+i • • • P''"-ipj„i(t) (4) 



1=0 
oven 



with P = i(P' + P), A = P' - P and f = A^. Here 
AS denotes symmetrization in i/, ...,/i„_i followed by 
anti-symmetrization in /i, v and subtraction of traces in 
all index pairs. The tensor operators are given by 



ASqia^^'iD^'^ ■■■iD^'^-^q (5) 



with D = {D — D)/2 and all fields taken at space-time 
point z 0. The analogous matrix elements of local 
vector quark operators are parameterized by ^^o(*) 
specified in For definiteness we consider in the fol- 
lowing A'!^'Q{t) and B^'^^^t) for up-quarks in a 7r+. Their 
counterparts for down-quarks and for tt~ or tt" readily 
follow from isospin invariance 0, since Wilson fermions 
preserve flavor symmetry. We note that A^^{t) is identi- 
cal to the electromagnetic pion form factor F-^{t), which 
we investigated in detail in . 

Lattice QCD results. — Based on our simulations with 
Wilson giuons and dynamical, non-perturbatively 0{a) 
improved Wilson fermions with Uf — 2, we have eval- 
uated the matrix elements in Eq. ^ for n = 1,2 and 
momentum transfers up to —t k, 2> GeV^. Configurations 
were generated at four different couplings /3 = 5.20, 5.25, 
5.29, 5.40 with up to five different k = Ksca values per 
13, on lattices of sizes F x T = 16^ x 32 and 24^ x 48. 
We have set the lattice scale a using a Sommer param- 
eter of ro = 0.467 fm Q. The pion masses are as low 
as 400 MeV, spatial volumes are as large as (2.1 fm)'^, 
and lattice spacings are below 0.1 fm (see Q for a list 
of lattice parameters) . The computationally demanding 
disconnected contributions present for even n are not in- 
cluded. For the tensor GFFs B^^g we expect them to be 
small in the physical limit, since they require a chirality 
fiip on a quark line and are thus suppressed by the quark 
mass jSj . All results were transformed to the MS scheme 
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FIG. 3: Study of discretization errors in _By'-|^(t=0)/m.^. 

at a scale of 4 GeV^ using non-perturbative renormal- 
ization 0. Further information on the computation of 
GFFs in lattice QCD can be found, e.g., in |6l. | l0l . and 
details of the present analysis will be given in [ll| . 

As an example we show in Fig. [2] the t dependence of 
-St(" =1,2)0 at /? = 5.29 and « 600 MeV. The ex- 
trapolation to the forward limit i = requires a param- 
eterization of the t dependence of the lattice results. As 
the statistics and t range of our data is not yet sufficient 
for sophisticated multi-parameter fits, we use a standard 
p-pole form F{t) — Po/[l ~ V(-P™p)]''i where the for- 
ward value Fq = F{t=Q) and the p-pole mass rup are free 
parameters for each GFF. Good fits are obtained in a 
wide range of p, with a preference for relatively low val- 
ues. On the other hand, a regular behavior of p^''{b±, s±) 
in the limit b± (which is of course inaccessible in 
a lattice calculation) requires p > 3/2 for B^'^Q{t) [1]. 
We therefore take p — 1.6 in the following. For the ex- 
amples in Fig. [2] we obtain BJ^'^^^{t—Q) — 0.856(60) with 
TOp = 0.949(57) GeV, and BJ^'^g{t^O) = 0.206(24) with 
TOp — 1.239(30) GeV. We stress that our final results 
show only a mild dependence on the chosen value of p. 
Taking, e.g., p = 2, which gives the power behavior for 
t — > — cx) expected from dimensional counting, changes 
our fits of Pj^o by less than the statistical errors even 



beyond the region —t<3 GeV^ where we have data 11| 
Before discussing potential discretization and finite 
size effects as well as the pion mass dependence of our 
results we note that, due to the prefactor in the pa- 
rameterization ffl, the GFFs B^^Q{t) must vanish like 
TOtt for TOtt ^ This is also required to ensure that 
the densities in Eq. ([T]) stay positive and finite in the chi- 
ral limit. In the following we therefore consider the ratio 
P^^jjg/mTT, which tends to a constant at = 0. 

Figure [3] shows the dependence of BJ^'^q (t—O) j niT^ on 
the lattice spacing a for two ranges of pion masses, where 
we have excluded those lattice data points which are most 
strongly affected by finite volume corrections (see below). 
We conclude that discretization errors are smaller than 
the statistical errors and neglect any dependence of the 
GFFs on a in the following analysis. 

Figure |4] shows the volume dependence of 
Pj'j^(i=0)/TO.n. for three different ranges of mT^. 
The finite volume corrections to the matrix elements 
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FIG. 4: Study of finite size effects of i?JjQ(t=0)/m.^. Shaded 
bands represent a combined fit (restricted to m^L > 3) in m.^ 
and L as described in the text. The dashed lines show the 
infinite- volume limit of the fit. 

with ri = 1 , 2 in Eq. ^ are known to leading order in 
chiral perturbation theory (ChPT) For m-^L 1 

the leading correction to i?J^Q(i=0)/m,r is proportional 
to TO^exp(— m^L) up to powers of {mT^L)~^^'^ , where L 
is the spatial extent of the lattice. Although our analysis 
includes pion masses as low as 400 MeV, we feel that a 
quantitative application of the chiral expansion requires 
lattice computations at even lower values of rnvr and 
probably the inclusion of higher-order terms. We take 
however the result of 12] as a guide to estimate the L 
dependence of our lattice data, fitting Bj'^'^ (t=0)/m.n- 
to the form cq + cim^ -I- C2m^ exp(— mTri). This 
fit, represented by shaded bands in Fig. |H gives 
By"o(i=0) = 1.47(18) GeV-i at L = oo and - 
440 MeV, compared to Bjfo(t=0) = 1.95(27) GeV^i 
at i ^ 1.65 fm as represented by the diamond in the 
lowest panel of Fig. |4l The typical corrections for 
B^'2Q{t=0) / rriTr are similar. Within present statis- 
tics, we do not see a clear volume dependence of the 
corresponding p-pole masses for n = 1,2. 

The pion mass dependence of B^'^Q(t=0) / 171^, is shown 
in Fig. m The darker shaded bands show fits based on 
the ansatz we just described. Data points and error 
bands have been shifted to L — 00. For — 140 MeV 
we obtain Bjfg (t=0)/TO^ = 1.54(24) GeV"^ with nip = 
0.756(95) GeV, and Bl^'^^{t=Q) / = 0.277(71) GeV"^ 
with TOp = 1.130(265) GeV, where in both cases we have 
set p — 1.6. The errors of the forward values include 
the uncertainties from finite volume effects. The light 
shaded bands in Fig. [5] show fits restricted to m^r < 650 
MeV using 1-loop ChPT [2] plus the volume depen- 
dent term C2m^ exp(— WTri). We note that the ChPT- 
extrapolation gives larger values for B^'^Q{t=0) at the 
physical point than the linear extrapolation in mj. 

To compute the lowest two moments of the density in 
Eq. ^ we further need the GFFs AIqII) with n = 1,2. 
For A^Q'{t) = F-^{t) we refer to our results in A 
detailed analysis of A^^{t) will be presented in [ll[, 
and first results are given in 
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FIG. 5: Pion mass dependence of B^'^g{t—0) / niTr ■ 
shaded bands represent fits as explained in the text. 
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p-pole parameterization with p = 1, which provides an 
excellent description of the lattice data and is consistent 
with power counting for t ^ —00. Fourier transform- 
ing the parameterizations of the momentum-space GFFs 
we obtain the densities p^{b±, s±). In Fig. [6] we show 
p"^^{b±, s±) for up-quarks in a 7r+ together with cor- 
responding profile plots for fixed bx- Compared to the 
unpolarized case on the left, the right-hand side of Fig. [6] 
shows strong distortions for transversely polarized quarks 
and thus a pronounced spin structure. The difference be- 
tween p = 1.6 and p = 2 for -BJj^q is negligible within er- 
rors. The negative values of the density on the lower right 
in Fig. [6l obtained for the maximal values of -Bj'j^ from 
the chiral extrapolations in Fig. [SJ are unphysical. They 
show that 1-loop ChPT cannot be regarded as quantita- 
tively reliable in this case and provides only a rough idea 
of the uncertainties related to the chiral extrapolation. 
From Eq. ([T]) wc obtain an average transverse shift 



J tPb^ b\p-^{b^,s^) 1 5?„o(t=0) 
/d25j_p"(6x,5±) " 2m, A-o(i=0) 



(6) 



We fit A^o"(0 to a 



in the y direction for a transverse quark spin s±_ = (1,0) 
in the x direction. Our lattice results give (6^)1 = 
0.151(24) fm and {b\)2 = 0.106(28) fm. 

Let us compare our results for B^^^ with those for 
the analogous GFFs BtuO that describe the dipole-like 
distortion in the density of transversely polarized quarks 
in an unpolarized nucleon. The corresponding average 
transverse shift is = BTno{i=^) / {2mN Ana{t—Q)) , 

where Ano{t—Q) is the n-th moment of the unpolarized 
quark distribution. With the lattice results of Q we 
find (6^)1 = 0.154(6) fm and (6^)2 = 0.101(8) fm for 
up-quarks in the proton. Remarkably, the distortion in 
the distribution of a transversely polarized up-quark is 
within errors of the same strength in a tt"^ and in the 
proton. An explanation of this finding has recently been 
proposed in the framework of quark models 13 1. 

The moments of the GPDs in the pion and Et in 
the nucleon can be connected with the respective Boer- 
Mulders functions, which describe the correlation be- 
tween transverse spin and intrinsic transverse momentum 
of quarks in an unpolarized hadron 1J|. They lead, e.g.. 
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FIG. 6: The lowest moment of the densities of unpolarized 
(left) and transversely polarized (right) up-quarks in a tt^ 
together with corresponding profile plots. The quark spin is 
oriented in the transverse plane as indicated by the arrow. 
The error bands in the profile plots show the uncertainties in 
^t'io (^~0) / rrin and the p-pole masses at m^^^'^ from a linear 
extrapolation. The dashed-dotted lines show the uncertainty 
from a ChPT-extrapolation (light shaded band in Fig. 5). 



the same sign and very similar magnitude as the cor- 
responding distortion in the nucleon Assuming the 
relation between impact parameter and transverse mo- 
mentum densities proposed in [3] this suggests that all 
Boer-Mulders functions for valence quarks may be alike, 
as argued in [l^. The large size of the effect might give 
new insight into the mechanism responsible for the large 
cos(20) azimuthal asymmetry observed in unpolarized 
np Drell-Yan production, which is sensitive to the prod- 
uct hi'^hi [i7|. It motivates future studies of azimuthal 
asymmetries in unpolarized irp and polarized np^ Drell- 
Yan production at COMPASS, the latter giving rise to a 
sm{(p + (j)s) asymmetry sensitive to hj^'^hi, where is 
the quark transversity distribution in the nucleon 

The numerical calculations have been performed on 
the Hitachi SR8000 at LRZ (Munich), apeNEXT and 
APEmiUe at NIC/DESY (Zeuthen) and BlueGene/Ls 
at NIC/FZJ (Jiilich), EPCC (Edinburgh) and KEK (by 
the Kanazawa group as part of the DIK research pro- 
gram). This work was supported by DFG (Forscher- 
gruppe Gitter-Hadronen-Phanomenologie and Emmy- 
Noether program), by HGF (contract No. VH-NG-004) 
and by EU I3HP (contract No. RII3-CT-2004-506078). 



to azimuthal asymmetries in semi-inclusive deep inelastic 
scattering (SIDIS) and in Drell-Yan lepton pair produc- 
tion. The density of quarks with transverse momentum 
k±_ and transverse spin in a 7r+ is given by 



f{x,k±,s±) = i 



(7) 



in terms of the unpolarized distribution /f and the Boer- 
Mulders function h^^. We notice the close similar- 
ity between ([7]) and the impact parameter density ([1]), 
but emphasize that k± and b± are not Fourier conju- 
gate variables. A dynamical relation between k± and 
b± dependent densities was proposed in 15] and implies 
h^'^ ^ —E^ for the distribution appearing in SIDIS — 
we recall that hj^"^ is time reversal odd and thus enters 
with opposite signs in SIDIS and Drell-Yan production 
With this relation, our results for -B|^„o imply that 
the Boer-Mulders function for up-quarks in a 7r+ is large 
and negative, and that its ratio to the unpolarized distri- 
bution is similar for up-quarks in a 7r+ and in a proton. 

Conclusions. — We have calculated the first two mo- 
ments of the quark tensor GPD in the pion. We 
find that the spatial distribution of quarks is strongly 
distorted if they are transversely polarized, revealing a 
non-trivial spin structure of the pion. The effect has 
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